function [a,len]=heap_sort(a,len)
% heap sort

[a,len]=heap_create(a,len);

n=len;
while(n>1)
    tmp=a(1);
    a(1)=a(n);
    a(n)=tmp;
    a=heap_adjust(a,1,n-1);
    n = n-1;
end
